Journal of Nonlinear Mathematical Physics 



Volume 10, Number 1 f2003l. IT^HsUl 



Article 



A Heat Transfer with a Source: 

the Complete Set of Invariant Difference Schemes 

Vladimir DORODNITSYN and Roman KOZLOV^ 

^ Keldysh Institute of Applied Mathematics, Miusskaya PI. ^, Moscow, 125047, Russia 
E-mail: dorod@spp.Keldysh.ru 

^ Department of Informatics, University of Oslo, 0371, Oslo, Norway 
E-mail: kozlov@ifi.uio.no 

Received February 28, 2002; Revised July 11, 2002; Accepted July 24, 2002 

Abstract 

In this letter we present the set of invariant difference equations and meshes which 
preserve the Lie group symmetries of the equation ut — {K{u)uj:)x + Q{u). All special 
cases of K{u) and Q{u) that extend the symmetry group admitted by the differential 
equation are considered. This paper completes the paper [J. Phys. A: Math. Gen. 
30, Nr. 23 (1997), 8139-8155], where a few invariant models for heat transfer equations 
were presented. 

1 Introduction 

Symmetries are fundamental features of the differential equations of mathematical physics. 
It yield a number of useful properties such as integrability of ODEs, symmetry reduction 
of PDEs, existence of various types invariant solutions, conservation laws for the invariant 
variational problems etc. Therefore, preserving symmetries in discrete schemes, we retain 
qualitative properties of the underlying differential equations. 

The purpose of this paper is to develop the entire set of invariant difference schemes 
for the heat transfer equation 

Mt = (K(n)u,)^ + Q(n), (1.1) 

for all special cases of the coefficients K{u), Q{u) which extend the symmetry group 
admitted by equation Hl.l|) . This paper is based on the Lie group classification j^I (see 
also [1]) of the equation (|1.1|) with arbitrary K{u) and Q{u). This classification contains 
the result of L V Ovsyannikov Q7] for equation (|1.1() with Q = as well as symmetries of 
the linear case K = 1, Q = 0, which were known by S Lie. 

A few examples of the invariant difference schemes and meshes were considered in [2] . 
In the present paper we complete the paper [2] going through all cases of K{u) and Q{u) 
identified in the group classification |5, we construct difference equations and meshes 
(lattices) which admit the same Lie groups of point transformations as their continuous 
limits. 
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Lie group analysis of difference equations is a very active field of research where many 
contributions were done, and various approaches were applied by several authors (see [18]). 
In our approach which we are following in this paper we pose the question: How does one 
discretize a differential equation while preserving all of its Lie point symmetries? Thus 
a differential equation and its Lie group symmetry are a priory given but not a difference 
model. One then looks for a difference scheme, i.e. a difference equation and a mesh, that 
have the same symmetry group and the same Lie algebra. The basic steps in this direction 
were done [5, 6, 7, 8, 9, 10, 11], which were summarized in a recent book The main 
idea is that the invariant difference equations and meshes can be constructed with the 
help of the entire set of difference invariants of the corresponding Lie group. In the next 
section we explain how to construct difference models that conserve the whole group of 
point transformations admitted by the differential equations. 

The article is organized as follows. Section 2 provides a brief overview of the invariant 
discretization procedure. In Sections 3, 4, 5 and 6 we consider the cases of an arbitrary 
heat transfer coefficient K{u), the exponential heat transfer coefficient e^, the power heat 
transfer coefficient u'^ and the special case of power heat transfer coefficient: u~'^^^ corres- 



Section 7 is devoted to the linear heat conductivity with a source. In particular, this 
section covers detailed study of the invariant difference scheme for the linear heat equation 
without a source {Q = 0) including such aspects as superposition principle, reduction of 
the invariant scheme on the optimal system of subalgebras and the way to transform the 
moving mesh scheme into a stationary one. Notice that in the paper [19] there were 
considered some difference approximations of the linear heat transfer equation, which 
preserve its different symmetries on different meshes. In [19] a difference equation and 
a mesh are a priori given, then it was shown that for some kind of the mesh there were 
preserved some symmetries of the linear heat equation and another meshes preserve other 
parts of the symmetries. Thus, there are no difference schemes which conserve the entire 
set of symmetries in one difference model. In Section 7 we will develop the difference mesh 
and difference equation, which conserve the complete set of original symmetries in the one 
and the same difference scheme. 

The same approach we will apply for some tens of other nonlinear models of heat 
equation (|1.1|) . Summarizing conclusions end up the consideration of the entire set of 
invariant schemes for the equation (|1.1() . 

2 Symmetry preserving discretization procedure 

1. Let us briefly describe this method, which was called the method of finite- difference 
invariants 0. 

Let the differential equation 



admit a known symmetry group G„, whose Lie algebra is spanned by the operators 
Xi , . . . , Xn of the form 



pondingly. 



F{t,x,u,ut,Ux, . . .) = 



(2.1) 




i = l,...,n. 



(2.2) 



18 



V Dorodnitsyn and R Kozlov 



where the coefficients ^* and rj are functions of t, x, u, since we consider Lie point 
symmetries. 

Then we would Hke to propose a discrete model 
F{z) = 0, 

n{z,h)=0, (2.3) 

where the first equation is the approximation of the initial differential equation and the 
second one defines a difference mesh. Both of these equations a priori are not given and 
we have to establish the invariant mesh, on which we should approximate the original heat 
equation. We will show, that in all special cases of equation (jl.lj) one can construct the 
system (|2.;-{|) . starting from the entire set of finite-difference invariants of the corresponding 
Lie group. 

We denote by z in (|2.3|) a finite number of difference variables, which are used in the 
considered difference stencil, i.e. a finite number of mesh points, which are needed for 
the approximation of the differential equation p.l|) . The equations (|2.3|) can be explicitly 
connected with each other (if invariant mesh depends on solution) or not. In the last case 
we can choose the invariant mesh firstly (for example, a fixed mesh) and then construct 
the invariant approximation of the original equation. If a mesh depends on the solution, 
all specifications of the mesh made in advance lead to restrictions on the symmetries which 
may be admitted by the considered discrete models. 

2. The idea of the method of finite-difference invariants springs from the invariant 
representation of differential equations. In the continuous case for the group Gn we can 
find the complete set of functionally independent differential invariants J = (Ji, J2, . . . , Jk) 
in the specified space which contains dependent and independent variables as well as the 
set of derivatives up to the highest derivatives involved in the formulation of the PDE ^5] . 
For the heat equation (|1.1|) we consider the space M ~ {t, x,u,ut,Ux,Uxx)- We prolong 
the operator ()2.2() on the variables of the space M 

ft fy (f 

pr X = X + C*TT^ + C"t7^ + 



with 



dut dux duxx ' 

C* = A(r/) - utDtie) - UxDtie), C = Dxiv) - utDxi^) - UxDxin, 

= Dxic) - utxDxie) - uxxDxin, 

where Dt and Dx are the total derivative operators for time and space correspondingly. 
Differential invariants are solutions of the system of linear equations 

pr Xi $(t, X, u, ut,Ux, Uxx) =0, i = 1,... ,n, 

and can be solved by standard procedure (see [15]). 

Then we represent the invariant differential equation in terms of these invariants 

F(Ji, J2,..., Jfc) =0. 

The obtained equation is invariant with respect to the group G„. 



A Heat Transfer with a Source: the Complete Set of Invariant Difference Schemes 



19 



3. In the discrete case the situation is more comphcated. Any given differential equa- 
tion can be approximated by means of infinitely many difference equations and meshes, 
which have the original differential equation as its continuous limit. The requirement of 
preservation of Lie group properties of the differential equation in its discrete counterpart 
still leaves some freedom in the approximations. Thus, as it can be seen recently, at some 
point we have to make a chose among the general family of invariant meshes. 

The structure of the admitted group essentially effects on construction of equations and 
meshes. Group transformations can break the geometric structure of the difference mesh 
that influences approximation and other properties of difference equations. First steps to 
the construction of the difference grids geometry based on the symmetries of the initial 
difference model were done in [5, 6, 7, 9]. There were found classes of transformations 
that conserve uniformity, orthogonality and other properties of the grids. 

It was shown [5, 6, 7, 9] that a transformation defined by 1)2. 2|) conserves uniformity of 
a mesh in t and x directions if and only if 



DDie) = 0, (2.4) 

-T — T 



DD{e)=0, (2.5) 

+h-h 

where D and D denote total difference derivatives in the time and space directions with 

±T ±h 

steps r and h correspondingly. 

For an orthogonal mesh to be conserved under the transformation, it is necessary and 
sufficient that 

Die) = -Din- (2.6) 

+h +T 

When condition (|2.Hj) is not satisfied for a given group, the flatness of the layer of a grid 
in some direction is rather important. For evolution equations it is significant to have flat 
time layers. There is a simple criterion of the invariance of flat time layers under the 
action of a given operator 1)2. 2|) : 

DDie) = 0. (2.7) 

±h+T 

These condition specify invariant geometry of grids for the given Lie group symmetries. 

If the operator coefficients do not depend on solution, then we can choose the 

invariant mesh as any solution of corresponding condition i|2.4|) ~ (|2.7j) . Otherwise the 
conditions (|2.4() - (|2.7() should hold on the solutions of the considered difference model. In 
that case we can figure a mesh out starting from the set of difference invariants. 

Further we choose a stencil which is sufficient to approximate all derivatives which 
appear in the equation. We will consider six-point stencils which have three points on 
each of two time layers. Such stencils allow us to write down both explicit and implicit 
difference schemes. For different transformation groups we will consider different meshes: 
orthogonal mesh which is uniform in space, orthogonal mesh which is nonuniform in space 
and nonorthogonal in time-space mesh, i.e. moving mesh. The corresponding stencils 
are different. Furthermore, the corresponding spaces of discrete variables are of different 
dimensions so that they have different number of difference invariants I = ili, I2, ■ ■ ■ , Ii) 
for the same Lie group G„. 
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For example, let us take an orthogonal mesh which is uniform in space. (We will 
describe later for which groups this mesh can be considered.) The stencil of this mesh is 
shown in Fig. 1. 

(x — /i, £, u_) {x,i,u) {x + h,i,u^) 



(x —li, t, 



U- 



ix,t,u) 



{x t,u+] 



Figure 1. The stencil of the orthogonal mesh. 

The corresponding discrete subspace is ten-dimensional: M ~ {t,x,T, h,u,U- ,u+,u, 
U-,u+), where t = t — t. The prolonged operator (|2.2() in this subspace has the form 



d 



d d 

+ V+ 



d 



d . d 
+ r] 



where we use time and space shifts notations / = f{t + t,x,u), /_ = f{t,x — h,u), 
= f{t, X + h,u). The number of functionally independent invariants is given by 



/ = dim M — rank Z, 



I > 0, 



(2i 



with dimM = 10 and the matrix Z composed by the coefficients of the prolonged on the 
space M operators 



z 



m 



yen itn-en) ' Vn 



Having found the finite-difference invariants as the solutions of system of linear equa- 
tions 

pr Xi ^{t, X, T, h, u, n_ , u+, U-,U-f-) =0, i = 1, . . . , n, 

we can use them to approximate the differential invariants 

Jj = fj{h,h, + 0(t", h"), j = l,...,k, 

where a and (3 define some fixed order of approximation. Notice, that approximation error 
0(r", h^) is invariant together with other terms in the above representation. Substitution 
of difference invariants li instead of differential ones Ji into the function F provides us with 
an invariant difference scheme. Practically we can often omit the representation of the 
differential equation in terms of its invariants and just approximate the original differential 
equation by the finite-difference invariants. The use of finite difference invariants is the 
main point in both ways. 
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So, the first step in the invariant approximation is the choice of the invariant mesh. 
The last step is the choice of the invariant discretization of the original equation on the 
invariant mesh. 

The described above method is algorithmic. We would like to stress that the invariant 
approximation in our way is still not unique. For example extending the stencil (means 
enlarging the number of mesh points involved in approximation) we can find invariant 
approximations of any higher order. 



3 An arbitrary heat transfer coefficient K{u) 

Now we start to develop invariant schemes going throw all cases of the Lie group classifi- 
cation [4]. Let us note that the group classification of the equation was done in |^ 
(see also jlj) up to equivalent transformations: 

62 _ g 

t = at + e, x = bx + f, u = cu + g, K = — K, Q = —Q, (3-1) 

a a 

where a, b, c, e, / and g are arbitrary constants, abc ^ 0. These transformations do not 
change differential structure of the equation (jl.lj) . transforming an admitted group into 
a similar group of point transformations. 

1. We start from general case, when the coefficients K{u) and Q{u) are arbitrary. 
Then the equation (|1.H) admits a two parameter group of translations only. This group is 
defined by the following infinitesimal operators: 

which generate the translations of independent variables. In this case almost no limits are 
imposed on a mesh and a difference equation. In particular, we can use an orthogonal 
grid in the plane {x, t) which is regular in both directions, as the conditions (|2.4)1 - (|2.6)1 
are valid for the operators 1)3.2(1 . 

The group with operators (|3.2() in the subspace {x,t,h,T,u,U-,u+,u,u- ,u-^-) corres- 
ponding to the stencil shown in Fig. 1 has eight invariants: 

r, h, u, n+, n_, u, n„, n+. 

That is why any difference approximation of the equation by the above invariants 

could give difference equation which admits the operators ()3.2() . For example, the explicit 
model 

where K{u) and Q{u) represent any approximation of the corresponding coefficients by 
invariants and Ux = '^x = right and left difference derivatives, admits the 

operators ()3.2|1 . 

2. If K{u) is arbitrary function and Q{u) = 0, the equation 



ut = {K{u)ux)x 



(3.4) 
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admits a three-parameter algebra of operators (see [T7]): 

This case is almost analogous to the previous one. The operators ()3.5() do not violate 
conditions of invariant orthogonality (|2.6)1 and invariant uniformity of a grid (|2.4|) . (|2.5j) . 
Thus, in this case we could use the orthogonal grid shown in Fig. 1. Any approximation 
of the equation (|3.4() by the seven invariants 

, U, U+, U-, U, U-, 

T 

gives an invariant model for the equation (jS3]). In particular the explicit scheme ^A.W^ 
with (5 = 0: 



r h\ \ 2 J h V2 J h ^ 

can be used. 

4 The exponential heat transfer coefficient K = 

In this paragraph we consider three cases of group classification for K = e^, in accordance 
with _4 and TT. 

1. If Q = then the equation 

ut = (e''n^.)x (4.1) 
admits a four-dimensional algebra of infinitesimal operators: 

d d d d d d 

As in the cases considered above, conditions of invariant uniformity and invariant or- 
thogonality are valid. A difference model for the equation ()4.1() can be constructed by 
approximation of the differential equation with the help of difference invariants: 

e"^, (n-u), (ii+-u), iu-u-), {u+-u), (u-U-). 

An example is the simple explicit difference model: 

exp — - — u a, - exp — - — u 5; , (4.3) 



r h \ ^ \ 2 J h 

but one has a lot of freedom to construct invariant schemes using the finite-difference 
invariants. 

2. For Q = 5 = ±1 we have a possibility to exclude the constant source from the 
equation 

Ut = (e"n^)^ + 5 (4.4) 
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by change of variables: 

u = u-6t, i = (5(e^*-l). (4.5) 

The equation H4.4() will be transformed into the equation H4.1() by this change, but the 
uniformity of the grid in the t-direction is destroyed. The equation (|4.4j) admits the four- 
dimensional algebra of infinitesimal operators: 

Xi = |, X, = ^, X, = e-''^+6e'''^, X, = x^ + 2^, (4.6) 
at ox at ou ox ou 

and we can easily see that the operator does not conserve uniformity of a grid in the 
time direction. The finite-difference invariants: 

-, {u — u — 6t), {u+ — u), (u — U-), (■U+ — n), [u — U-) 



permit us to construct the following variant of difference model for the equation ()4.4() : 
6(u-u)-T 1/ /u+ + u\ fu + u^\ \ 

Let us note that the change of variables ()4.5I) transforms the model ()4.7|) considered on 
the orthogonal grid with the time interval [0,T], given by the formula 

t„ = 51n(l + ^(e^^-l)) , n = 0,...,k, (4.8) 

where k is the number of time steps of the grid, into the model (|4.3|) with uniform time 
grid on the time interval [O, 6 (e*^^ — l) ] . 
3. If Q = ±e"", a 7^ the equation 

ut = (e«n.). ± e"" (4.9) 

admits 3 infinitesimal operators: 

These operators satisfy the conditions of orthogonality and uniformity of invariant grids 
and we will consider the stencil of Fig. 1. Any approximation of the equation 1)4. 9|) by the 
finite-difference invariants 

a-l 
T 2a 

— ; — , e°"r, (u — u), (ii+— n), (u — u-), (u+ — n), (n — n_) 
h 

gives a variant of a difference model for the equation (|4.9|) . admitting the symmetries 
(|4.inj) . for example, we obtain the following model: 
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4. In accordance with group classification we will also consider the case Q = ibe"+5, 
(5 = ±1. As in the case 2 we have the possibility to exclude the constant source from the 
equation 

= (e"n,.)x ± + (5 (4.12) 

by the change of variables (|4.5|) . The equation (|4.12|) will be transformed into the equa- 
tion ()4.9|) . The equation ()4.12|) admits the following infinitesimal operators: 

Finite-difference invariants of (|4.13j) 

e^{e^'^ — l), h, {u — u — 6t), (n+ — u), (u — U-), {u+ — u), {u — ii-) 
permit to construct the following variant of the difference model 

- 1 (- (^) > - - (^) ^) 

on the invariant grid (|4.8|1 . The model for the considered equation can be obtained from 
the model (|4.1H) with the help of the transformation (|4.5|) . 

5 The power heat transfer coefficient: K = u^, cr 7^ 0, — | 

For K = u"^ further classification depends on the source. 
1. Let us start with the simplest case Q = 0: 

Ut = {u''Ux)x- (5.1) 

Symmetries of the equation ()5.1|) are described by the four-dimensional algebra of 
infinitesimal operators (see jl7j): 

d d d d d d , , 

^1 = 7^, ^2 = 7r> ^3 = 'it— + x—, Xi = ax— + 2u—. 5.2 
ot ox at ox Ox ou 

With any a the operators ()5.2|) conserve uniformity and orthogonality of a grid. The 
finite-difference invariants corresponding to the stencil of Fig. 1 are 

^ r u n_ n_ 

a"^ u u u u u 

They permit us to write, for example, the following variant of the difference model on the 
orthogonal uniform mesh in both directions: 

Orthogonal grid is not the only possible way for discrete modeling. On the example 
of the equation 1)5. 1|) we will show how to introduce a moving mesh of the form shown in 
Fig. 2. 
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t 




Figure 2. A moving mesh with flat time-layers. 



One can use an adaptive grid defined by an evolution equation (see also jS]) 



dx 

— = (p{t,X,U,Ux). 



(5.4) 



In this case the heat transfer equation will take the form 



dt 



(5.5) 



Different requirements could be imposed on the function (p. If we require invariance 
of the equation 1)5. 4|) with respect to the whole set of the operators (|5.2|) , our freedom to 
chose If is limited by the function 



C = const. 



Below we show how to introduce Lagrangian type of evolution ^ . Let us note that the 
equation ()5.1|) has a form of the conservation law that presents the conservation of heat. 
Hence we can search for a moving mesh of Lagrange type which evolves in accordance 
with heat diffusion. We should find an evolution ^ which satisfies the equation 



d r'^^^^ 



dt 



udx = 0. 



xl{t) 



Since 



dt 



x2 



udx 



xl 



x2 



xl 



du 
'di 



dx + 



dx 
I— 
dt 



-I x2 



xl 



U Ux + U 



dx 
'dt 



x2 



xl 



we obtain the evolution ^ 



,<T-1 



Ux- Note that this evolution is invariant with respect 



to the operators (|5.2j) . Our initial differential equation l|5.1j) can now be presented in the 
form of the system 



dx 

— = -U Ux, 

dt 



du _ , ^, , 

— = -u uxx + {(y - 



u^. 



(5.6) 
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Let us mention that the equation (|5.1|) has two conservation laws 
Ut = {u'^Ux)x and {xu)t = I xu'^Ux 



a + 1 

For the evolution system (|5.6() it is convenient to present the conservation laws in the 
integral form 



u dx = 0, — / xudx 



dt Jxi(t) dt J^m^ cr + 1 



X2 



(5.7) 

For difference modeling of the system (|5.6() we can consider the stencil shown in Fig. 3. 



{x — h ,t,u^) {x,t,u) (x + /i+,t,n+) 

{x — h~,t,u-) {x,t,u) {x + h'^,t,u+) 

Figure 3. The stencil of the evolutionary mesh. 

In the space of the variables {t, x, r, h'^ , ,h'^ , , Ax, u, Uj^,U-., u, u+,U-) correspon- 
ding to this stencil there are ten finite-difference invariants: 

r u u+ U- u_ Ax 

u 



Approximating the system ()5.6() by these invariants we can get, for example, the system 
of two equations: 

where we approximated the heat conservation law to obtain the equation for the solution u. 

The first equation of system (|5.6|) shows that the evolution of x depends on the solu- 
tion. The system (j5.8|) may be inconvenient for computations because steplength will be 
changed automatically and the nature of this process is not clear. In order to avoid this 
indeterminacy we introduce a new space variable which values characterize the evolution 
trajectories of x. Let us consider the variable s defined by the system: 

St = U^Ux, Sx = u. 

It is easy to see that each trajectory of x is prescribed by a fixed value of s since 
ds dx ^ ^ , n\ r. 

-TT = St + Sx-JT =U Ux- U-{U )x = 0. 

dt dt a 

In a new coordinate system with the independent variables (t, s) the equation ()5.1() has 
the form 
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and the former space variable x satisfies 

Xt = —u'^Usi Xs = — . (5.10) 

u 

For discrete modeling of the equation (|5.H) one can use the equation (|5.9j) in the new 
independent variables (i, s) to describe the diffusion process and the first equation of the 
system (|5.1()|) to trace the evolution of the coordinate x. The equation (|5.9|) considered 
together with the system (|5.10p admits the following symmetries 

d d d ^ d d d 

at ox OS ot OS ox 

X5 = (. + 2)4 + .xA + 2.|-. (5.11) 

In the new variables (t, s) the stencil becomes orthogonal so that there is no need to 
consider a nonuniform grid in the variable s. There are following invariants for this set of 
operators in the space (t, t, s, hs,x,h^ , h~ , , h~ , Ax, u, u^,u^,u, «-!-,«_) corresponding 
to the orthogonal stencil in (i, s) extended by additional dependent variable x: 

^ T u u+ u+ U- h~ h~ /i+ Ax hg 



ht"^ ' u M ' M ' u ' ii ' hf hf hf ht ' ht 

By means of these invariants we get an approximation of (|5.9|1 which has the form of 
a conservation law 

-(--- 




(5.12) 

where < a < 1. Note that in the coordinates (t, s) variable x is introduced by the 
system (|5.10p as a potential for the equation (|5.9j) . Similarly we can introduce x as 
a discrete potential with the help of the system 

Ax a (ul+^-u'L+^\ 1-a 



r o- + 1 \ 2hs a + l \ 2hs I 

hs 2 \u u+ J 

In computations only the equation ()5.12() and the first equation of H5.13() are needed. 
The second equation of the system (|5.13|) is needed only to establish the connection bet- 
ween solutions u{x) and u{s) for a fixed time. Given some initial data u{0,x) = uq{x), we 
choose an appropriate steplength hg for the Lagrangian mass coordinate s. Then we can 
introduce the mesh points Xi in the original coordinates satisfying 

^i+l 1 I ^ 



2 \uo{xi) uo{xi+i] 
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i.e., we use this equation to establish difference relation between the original space coor- 
dinate X and the Lagrangian mass coordinate s. Computing the solution with the help of 
the numerical scheme (|5.12|) and the first equation of we preserve the relation 

Xj+i - Xi _ 1 f 1 ^ 1 



/is 2 \ Uj u-i+i 



Introducing the mass type variable s, we can rewrite the conservation laws 1)5. 7p as 
as = 0, — / xas 



dt Jgi dt Jgi a + 1 

The proposed discrete model possesses difference analogs of these conservation laws 



hg = const, 



i=l 

N-1 . , . Af-1 , I „<^+l\ 

E Xj + Xj+i , _ Xj + Xj+i _ a I ^AT+l + \ 

2 '2^ 2 '~cr + ll 2 ] 

i=l i=l \ / 

1-a ( uU,\+uV-' \ a l-g / Cr+^o+^ \ 

a+iy 2 J a+iy 2 J a+l^ 2 J ' 

Let us mention that for computations we need to propose some method for the boundary 
points. 

2. Q = 5u, 6 = ±1. In this case the symmetry of the equation 

ut = {u''ux)x + 5u (5-14) 

is described by the following infinitesimal operators: 

^ d d ^ d 

Xi = —, X2 = — ' ^3 = crx— + 2n— , 
at OX ox Ou 

= (5.15) 

By the change of variables 

u = ue~^\ t = -(e'^"*-l) (5.16) 
a 

the equation (|5.14|) is transformed into the equation (|5.H1 . 
The finite-difference invariants 

(e^"'' - 1) u \ Uj, u- u+ u. 
In r 



' \ u y' n' u' -u' ii 
give the following possibility for an explicit difference model 

^^ln^-rl=iff^Vn.,-f^) n,). (5.17) 



S'^-r - l\ u J h\\ 2 h'' \ 2 h 
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Let us remark that the change of variables transforms this equation considered on 

the orthogonal mesh with time layers 

t„ = ^ln(l + ^(e^'^^-l)) , n = 0,...,k, (5.18) 

which fill the time interval [0,T], into the equation ()5.3j) on the uniform grid on the time 
interval [O, |(e'''^^ - l)]. 

3. Q = ib'u'^+i + 5n", 6 = ±1. The equation 

ut = {vFux)x ^ u"" , cr,n = const, (5.19) 

admits a three-parameter symmetry group. A possible representation of this group is by 
the following infinitesimal operators: 

The set 1)5. 20() satisfies all the conditions (|2.4|) - H2.6|) . So we can use an orthogonal grid 
that is uniform in both t and x directions. By considering the set of the operators (|5.2U|) 
in the space {t,t, x, h'^ , h~ ,u,u^,U-,u,u^,U-) that corresponds to the stencil shown in 
Fig. 1 we find 7 difference invariants of the Lie algebra: 



t2("~i) u n_|_ n_ n+ U- 

n u u u u u 

The small number of symmetry operators provides us with a large number of difference 
invariants. Thus we are left with some additional degrees of freedom in invariant difference 
modeling of 1)5. 19(1 . By means of the discrete invariants we obtain the following explicit 
scheme: 

where n~ = "+.7" , u 



h ^ h ^ 
4. Q = ±u'^^^ + 5u, 6 = ±1. The equation 

ut = {u''ux)x ± u^'^^ + 6u (5.22) 

is connected with the equation ()5.19|) by the transformation ()5.16|) . The infinitesimal 
operators admitted by the equation are 

With the help of invariants for the operators 1)5. 23(1 

u [e -1), h, [5\n r , — , — , — , — 

\ u J u u u u 

we have the following example of an explicit difference model: 

5\n--T] = - [[^ — Ux-\ — - — Us ±n'^+\ (5.24) 



_i\ u / h\\ 2 h \ 2 h 



This equation considered on the grid (|5.18|) is transformed into the equation 1)5. 21() on 
a uniform time grid by the variable change 1)5. 16() . 
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6 The special case of power heat transfer 
coefficient: K = u~^^^ 

1. If Q = 0, then the symmetry of the equation 



ut 



-4/3 



U 



(6.1) 



is described by the five-dimensional algebra of infinitesimal operators (see |17j): 



5t ' 5x ' 

2x- 3'u— , Xs 

ox ou 



2t— X— 

dt dx' 



d 



d 



X — 3xu-—. 



dx 



du 



(6.2) 



These operators conserve orthogonality and uniformity of a grid in the time direction. The 
operator conserve uniformity in the i-direction, but does not conserve uniformity of 
the grid in the j;-direction; however orthogonality is not disturbed. We will consider the 
stencil shown in Fig. 4. 



{x — h ,t,u^) {x, t, u) 



{x + h+,i, u^] 



{x — h ,t,U-) {x,t,u) (x + /i"*", t, u+) 

Figure 4. The stencil of nonuniform mesh. 
The finite-difference invariants corresponding to this stencil 

u u+ n_ 1/3 h+ 1/3 ^3 h- v?l'^ ( h+h- 
u u+ u- \Jt ^Jt vr \n+ + n 

^ive among others the explicit difference model 

u-u ^ h+ + h- f u-^'^-u-^/^ _ n-V3_^:V3 \ 
r ~ Qh+h- \ h+ h- ]' 



(6.3) 



Remark. Let us note that we can not propose a space mesh /i+ = f{x,h-) which is 
preserved under all transformations of the group ()6.2|) . It can be clearly seen from the 
absence of difference invariants in the space (x, /i_ , /i+). For example, if we take a solution 
of the difference scheme 1)6. 3() on a regular mesh h- = /i+ , a general group transformation 
corresponding to ()6.2() will transform the solution into another solution of this difference 
scheme but possibly on a nonuniform mesh. This remark is also valid for the cases 2, 4 
and 5 of this section. 

2. In the case Q = 6u, 5 = ±1 equation 

ut = {u-''/^u^)^ + 6u (6.4) 
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admits operators 

d d d d 

^1 = ^' ^2 = ^3 = 2x- 3u— , 

at ox ox ou 

ist d ist d odd , , 

X4 = e 3 — + 5e 3 n— , X5 = x^—-3xu—. (6.5 
ot ou ox ou 

By the change of variables (|5.16|) this equation can be transformed into the equa- 
tion (|6.H) . Let us write out the difference invariants for the set of the operators (|6.5|1 : 

ii \ 2/3 ^ h'^h" \ 1 
5 In T , u^''^ 



u 



h+ + h-J ^{e&^T _ Y) 



u 



^J\y^h+ u'J\y^h- u]/\^/^h+ u^'\y^h- 



^/{e^^^ - 1) ' Vle"^"" - 1) ' ^(e-*"" - 1) ' 

These invariants can be used for construction of a difference model for the equation (|6.4() . 
We show the explicit variant of the difference model: 

au /^u \ h+ + h- (u-'/'-u-^/^ n-i/3_^:i/3\ 

'^i^r-^ = - TT rz • (6-6) 



eS<^T _i y- J Qh+h- \ h+ 

3. Q = zbti", n 7^ — |. The equation 

ut = {u-^'\,)yu^ (6.7) 
admits infinitesimal operators: 

Although this equation is specified in group classification (see j4j), it is a particular case 
of the equation ()5.19p . — there is no extension of the admitted group. That's why as 
an invariant difference model for the equation (|6.7() we can use the model ()5.21() with 
parameter o" = — |, corresponding to the given equation. 

4. If Q = au~^/^ , a = ±1, then the variant of the difference model for the equation 

= (n-4/V),±n-i/3 (6.9) 

depends on the sign of the coefficient a. The equation 1)6. 9() admits a five-dimensional 
algebra of infinitesimal operators, namely 

d ^ d A d ^ d 

Xi = — , X2 = —, X3 = -t— + 2n— , 
ot Ox 3 ot Ou 

ox Ou 

X5 = e-2\/^-^ + V3^e-2V^-^i_. (6.10) 
Ox Ou 

a.) The case a = 1. By the change of variables 

u = ucosh^ ( —pz I , X = \/3tanh ( —pz ] (6.11) 
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we transfer the considered equation into the equation (|6.H) (see ^). With the help of 
difference invariants 

U UA- ?i_ r,,n { 1 1 



« \tanh(^j 



-^M^/^uV^ sinh ( — T= \ , ^^u^^^u^^ s'mh 



h- 

71 



we can construct a difference model. Let us write out one of the possible variants of the 
difference model, namely an explicit model: 

u — u 1 / 1 1 

+ 



tanh(^) tanh(^) 



sinh (^) sinh 



(6.12) 



b.) The case a = — 1. By the change of variables 



u = u cos 



' (^) ■ ^ = (^) <"^' 



we can transfer the given equation into the equation (|6.1jl (see [J)- The set of finite- 
difference invariants: 

u n_i n_ ^ r,/o / 1 1 
/rn"^/^ 7 ^ + 



^ ^- Itan(^) tanf^ 



,\/3y Vv^/ 

1 1/3 1/3 . 1 1/3 1/3 . fh-\ 

u ' u_i sm —= , —^u ' sm 



- 



provides us with a possibility to construct an invariant difference scheme. For example, 
one can use an explicit difference model: 

u — u 1 / 1 1 

+ 



:V3 _ ^-1/3 / Ai) ^-1/3 COS ('^ j _ ^:i/3 ^ 



(6.14) 



We stress that the obtained difference models (|6.12|) and (|6.14j) are connected with the 
difference model 1)6. 3() for the equation 1)6. 1|) by the changes of variables ()6.11() and (|6.13p 
correspondingly as the initial differential equations. 

5. Q = au~^/'^ + 5u^ \a\ = \5\ = 1. As in the previous point, two cases of parameter a 
in the equation 

ut = {u^'u^)^ ± n''+^ + 5u (6.15) 
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should be considered separately and two difference models should be constructed. Let us 
write out the infinitesimal operators, admitted by the equation (|6.15|1 : 

d d ist d ^ 4M d 

^1 = 777, X2 = —, = e 3 — + <5e 3 u—, 

at ox at ou 

ox ou 

= e-2^^-|- + Vs^e-'V^'-u^. (6.16) 
ox ou 

a.) The case a = 1. The change of variables 1)5. 16() transforms the considered equation 
into the equation 1)6. 9|) and the change ()6.11|1 into the equation H6.1() . 

We write out the set of finite-difference invariants for the equation (|6.15j) with a = 1: 



u ) Vtanh(^) tanh(^^ 

u^^u^f /h+\ u^/^uj^ . fh~\ 

— — sinh —= , — — sinh —= , 

u'/^ii'f , fh+\ u'/^u]/^ . fh~\ 

sinh —;= , — — sinh — — 



^^^sar _ 1) VVsy ' - 1) v\/3y ■ 

The explicit variant of the difference model for the equation H6.15() on the time grid (|5.18l) 
has the form: 

au f ^, u \ 11 1 



/ u;'/»-„-/^cosh(^) ,.-i/»cosh(^)-C'^° \ ^^^^^ 

b.) The case a = —1 Using the change of variables 1)6. 13() we can transfer this equation 
into the equation (|6.4j) and by the change (|5.16j) into the equation (|6.9|) . Difference model 
for the equation (|6.15|1 can be obtained with the help of the invariants 



(t^) tan 



sm —;=z , — = sm 



V(e'5'^^ - 1) V\/3/' - 1) V\/3/' 



: sm I — = . — — sm , _ 



34 



V Dorodnitsyn and R Kozlov 



One of possible difference models for the equation (|().15j) is 




(6.18) 



The difference models ()6.17|) and ()6.18|) are connected with the model ()5.6() by changes 
of variables ()6.11|) and H6.13() correspondingly. The variable change 1)5. 16() transforms 
the difference models obtained in this point into model of the point 4 for corresponding 
values of the parameter a. This example shows that in invariant difference modeling it is 
possible to get consistent models which are connected with each other by the same point 
transformations as their original differential counterparts. 



7 Linear heat conductivity with a source 

In this section we consider the semilinear heat transfer equation 

Ut = Uxx + Q{u) (7.1) 

with different types of a source. 

1. With Q = zize" the equation becomes 

ut = u^x ± e". (7.2) 
It admits a three-dimensional algebra of infinitesimal operators 

d d d d d 

^1 = 7^, ^2 = 7^, X^ = 2t- + x—-2—. (7.3) 
at ox ot ox Ou 

It is easy to check that for the operators the conditions of orthogonality and uniformity 
conservation of a grid hold. Approximation of the equation by the invariants 

— , re", {u — u), {u+ — u), (n — u_), (u^ — u), (u — u_) 
r 

will give different types of difference models. An explicit one is 

^^ = Uu^-Us)±e\ (7.4) 

T tl h h 

2. Q = ibu". The equation 

Ut = Uxx ± u"- (7.5) 
admits the following infinitesimal operators: 
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The operators satisfy to conditions (|2.4j) - (|2.()j) and the finite-difference invariants 

/i^ u Uj^ ti_ M4. n_ 

T U U U U U 

permit us construct, for example, the following difference scheme: 

'^ = Uu^-u^)±u-. (7.7) 

T n h h 

3. Q = 5tilnn, (5 = ibl. The semilinear heat transfer equation 

ut = Uxx -\- Sulnu, 5 = ±1 (7.8) 

admits the four-parameter Lie symmetry group of point transformations 0] corresponding 
to the following set of infinitesimal operators: 

Xi = |^, ^2 = 1-, X^ = 2e''^-6e''xu^, x, = e''u^. (7.9) 
ot ox ox ou ou 

Before constructing a difference equation and a grid that approximate (|7.8j) and inherit 
the whole Lie algebra 1)7. 9|) we should first check condition ()2.6|) for the invariance of 
orthogonality. The operators Xi, X2 and X4 conserve orthogonality, while X3 does not: 
the condition 1)2. 6|) is not true for operator X3. Consequently an orthogonal mesh can not 
be used for the invariant modeling of (|7.8j) . Conditions l|2.7|) are true for the complete set 
of operators, so it is possible to use a nonorthogonal grid with flat time layers and we will 
use the grid shown in Fig. 2. 

A possible reformulation of equation 1)7. 8() by using the four differential invariants in 
the subspace {t,x,u,Ux,Uxx,dt,dx,du): 

/ "Ux \ ^ Uxx T Ux dx 

Ji=dt, J2= J3 = 2^ + — , 

\ u / u u at 

du 1 f dx\^ 

J 4 = — ; in n + -- — 

udt A\dt J 

is given by the system 

J3 = 0) J'i = J2 
that is 

dx du ^ , ui 

— = -2—, — =Uxx + Sulnu-2^. (7.10) 
dt u dt u 

So, the structure of the admitted group suggests the use of two evolution equations. 

As the next step, we will find difference invariants for the set X1-X4 of the group ()7.9() . 
These invariants are necessary for the approximation of system ()7.1U() . We will use the 
six-point difference stencil of Fig. 3 on which we will approximate system ()7.1U() . The 
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stencil defines the difference subspace {t,t,x,x,h^,h ,h^,h ,u,uj^,u^,u,uj^,u-) and 
the group (|7.9|) has the fohowing difference invariants 

h = T, l2 = h^, I^ = h', I^ = h^, 1^ = %", 
Iq = (Inu)a; - (Inn)j;, I7 = (In'u)a; - (Inu)j;, 

I9 = 6 Ax + 2(1 - e^^^) I . ^ . (hin)^ + . (Inn)^ 

y /i+ + /i- h+ + h- 

2 , /in „-5tn A„ J 



ho = 5{AxY + 4(1 - e-'^) {\nu- e'^ In-uj , 

1 A ^ /IN \nu+—\au r-i \ Inn— lnii_ 

where Ax = x — x, (Innj^; = — ^ , (inuj^ = . 

An explicit model can be chosen 

-'8 — U, iio — 7 — J — —f i6, 

h + h 

i.e. 

i(Ai)' + 4(1 - e-^) (ini - e'^ Inu) = ? *f!lli^ [(lnu)„ - (Inti)j]. (7.11) 

One invariant solution of this scheme is given in [2] . 

4. A linear heat equation without a source (Q = 0). 

4.1. Prehminary consideration. The linear heat transfer equation 

ut = Uxx (7-12) 

admits a six-parameter Lie symmetry group of point transformations, corresponding to 
the infinitesimal operators 

d d d d d d 

Xi = Tr. ^2 = 7^, ^3 = 2t^^ xu— , X4 = 2t— +x— , 

ot ox ox Ou ot ox 

and an infinite-dimensional symmetry 
X*=a{x,t)-^, 

where a{t,x) in an arbitrary solution of equation 1)7.12(1 . Symmetry X* represents linearity 
of the equation ()7.12|) . 

Probably, the simplest approximation of the linear equation is the explicit scheme 

u — u u+ — 2u + U- , ,N 
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considered on a uniform orthogonal mesh. This equation is invariant with respect to the 
operators Xi, X2, and of the set (|7.13() . Since the equation is hnear it possesses 
a superposition principle that is reflected in the invariance with respect to the operator 

d 

XI = ah{x,t) — , 

where ah{x,t) is an arbitrary solution of equation 1)7. 14() . In [2] it was shown how to 
construct a discrete model which admits the six-dimensional group (|7.13|) . 

To preserve the Galilean operator X^ and the projective operator X^ it is necessary to 
introduce a moving mesh. 

4.2. Heat transfer system of equations and superposition principle. With the 
help of the differential invariants of the operators 1)7. 13(1 in the space {t,x,u,Ux,Uxx,di, 
dx, du) 

dx + 2^dt ^ du Idx^ ( Uxx ul\ , 

we can represent the heat equation (|7.12|1 as the system 

Jl = 0, J2 = 
that is 

dx Ur du ul , , 

— = -2—, —=Uxx-2—. (7.15 
dt u dt u 

By construction this system is invariant with respect to the six-dimensional group gene- 
rated by the operators (|7.13|) . The system also inherits the superposition principle of the 
linear heat transfer equation. The superposition principle has the form of summing two 
solutions of the system ()7.12|) . but it also acts on the trajectories on which the variable x 
evolves. For two arbitrary solutions Ui{t,x) with the trajectories Xi{t): 

dXi ^ 2Uix 
dt Ui ' 

and U2{t,x) with the trajectories X2{t): 

dX2 _ 2U2X 

dt U2 ' 

their linear combination 

U = aUi + PU2, a, 13 = const, (7.16) 

is also the solution of ()7.15|) . However, this linear combination has its own trajectories 
satisfying 

dX _ aUi dXi PU2 dX2 

dt aUi + (3U2 dt aUi + (3U2 dt ' ^ ' ' 
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Therefore, the superposition principle can be presented in the following form: 




(7.18) 



Let us show the superposition principle for system (|7.15|1 by an example. The solution 

1 ( {x-af 

Ui = —== exp ' 



has the trajectories 

.t + h 

X = a + [xo - a)— — , 
while the solution 

U2 = . ^ exp ' ' 



exists on the trajectories 

t + t2 



X = b + {xq — b)- 



t2 



The linear combination ()7.16() of these two solutions is also the solution of system 1)7. 15() . 
Its trajectories are 

dX _ ^/t+tl \ 4(t+ti) ) \t+ti ) ^ ^Tt+t^ I, 4{t+t2) J \t+t2 



dt " cxpf I ^ cxpf 

Examples of evolution of grid points and corresponding solutions are shown on Figs. 5-8 
(for computations we used discrete model ()7.19|) which will be introduced in point 4.3 of 
this section). 

If we consider Lagrangian derivatives of the solution f7, we could get the superposition 
principle 

du ^dx\ (dUi dxA , /dU2 dX2 

— - U:r.— I = a I —r. Uix—^ ] + P \ —rr - U2x 



dt dt J \ dt dt J \ dt dt 

i.e. the superposition principle Ut = aUu + (3U2t of the linear heat equation expressed in 
terms of the of total derivatives of U and X. 

4.3. Invariant schemes on moving meshes. For the difference modeling of sys- 
tem ()7.15() we need the whole set of difference invariants of Lie symmetry group ()7.13|) in 
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-20 -15 -10 



10 15 20 



Figure 5. Evolution of solution (|7.16() : a = /3 = 1, ii = ^2 = 10, 



5 = 8. 





Figure 6. Mesh trajectories for the solution shown in Fig. 5. 
the difference space corresponding to the chosen stencil (t, t, x, x, , , h'^ , h~ , u, it, u+, 



h 
h 
h 
h 
Is 



h+ ^ h+ ^ h+h+ 
h h- r 



T^l'^ii /l(Ax)^ 

— — exp 

/i+ M V 4 r 



1 /^^ 



4 r /i+ + /i V /i+ u /i~ u 



1 /i+2 
+ 



1 , 1 , n_\ 

— In^ + ^In— , 



4 r h+ + h- \h+ u h- u 

Axh+ 2h+ fh~ u+ h+ 

+ — — I —In— - — In- 



Axh+ 2h^ 
+ 



- — In— — In— — 



T h+ + \h+ u h" u 
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-20 -15 -10 



10 15 20 



Figure 7. Evolution of solution lf7?TB|l : a = 0.25, /3 = 1, ii = = 10, a = -8, h 





Figure 8. Mesh trajectories for the solution shown in Fig. 7. 



Approximating system (|7.15j) by the invariants, we obtain a system of difference evo- 
lution equations. As an example, we present here an invariant difference model which has 
explicit equations for the solution u and the trajectory of x: 



Ax 



2r 



h+ + h- 
2 / 1(A2;)2 
2 r 



h+ u u 



-] exp 



1 



4t 



We also can write out an implicit model 



h+ u a u 



(7.19) 



Ax 



2r 



h , u+ ^ 
-—In— + —In — 

h+ u h~ u , 



1 + 



4r 



1 ti_ 1 u_ 

— In h — In — 

h+ u h- u 
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It is also possible to combine an explicit equation for the mesh and an implicit approxi- 
mation of the PDE or vice versa. Other ways to approximate the system 1)7. 15|) are also 
possible. 

4.4. Optimal system of subalgebras and reduced systems. Among all invariant 
solutions there is a minimal set of such solutions, called the optimal system of invariant 
solutions. From this set of invariant solutions any invariant solution can be obtained 
by an appropriate group transformation. The difference model (|7.19|) is a system of two 
evolution equations. To find its invariant solutions we need to provide a time mesh which 
is invariant with respect to the considered operator. An invariant time mesh giving flat 
time layers can be represented by the equation 

n = 9{ti), i = 0,l,2,.... (7.20) 

We request this equation to be invariant with respect to the considered symmetry. Since 
for the operators 1)7. 13(1 the coefficients i^* do not depend on x and u we can propose an 
invariant time mesh for any symmetry. In the case ^* = the function g can be taken 
arbitrary. For example, we can choose the uniform mesh tj = jr, t = const. Thus, 
different invariant solutions may have different time meshes. 

The adjoint action of the Lie group transforms an invariant solution into another 
one |15| I14j . In our case it also transforms the time mesh equation (|7.2fl|) . Thus the 
adjoint action gives us a new invariant solution with a corresponding invariant mesh. 

On the example of the difference model (|7.19|1 we will construct the optimal system of 
solutions which are invariant with respect to one-parameter groups. The optimal system 
of one-dimensional subalgebras of the algebra of symmetries for the linear heat equation 
consists of algebras corresponding to the operators (see [Tl|): 

OX ou ot ou 

d d d d d d 

Yi = Xi-X^ = — - 2t— + xu—, Y5 = X4 + 2cXe = 2t— + x— + 2cu—, 
ot ox Ou ot ox Ou 

^6 = ^1 + ^5 + cXe = (4t2 + 1)_ + 4tx— + {c-x^- '^t)u—. 

Let us find invariant solutions corresponding to these one-dimensional subalgebras. 

1) The subalgebra corresponding to the operator Yi has only constant solutions u = C, 
C = const considered on an orthogonal mesh Ax = 0. 

2) The subalgebra corresponding to the operator Y2 does not have invariant solutions 
(the necessary condition of existence of invariant solutions does not hold llSj). 

3) The operator 13 has the following invariants: tiexp(— ct), r and Ax. The time 
step T is invariant, so we can consider a uniform time mesh. We will seek a solution of the 
difference model in the form 



u = exp(ct)/(x). 
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Substituting this invariant form of the solution into system (|7.19|) . we get: 



h+ + h- \ h- V fix) J h- \ f{x 

fix) V / lAx' 
' exp — 2cr 



f{x + Ax)J V 2 r 

System 1)7. 2ip will become a system of two ordinary difference equations if we project 
it into the invariants space. To project the system we have to impose 

Ax = -h', 0, or (7.22) 



A solution of system (|7.2H) with one of conditions (|7.22|) provides the solution of sys- 
tem (|7.19p which is invariant with respect to the operator Y^. 

4) The operator I4 has invariants: uexp (— — ^t^), x + t"^, r and ^ — t. Let us 
search a solution of the difference model ()7.19|) in the form 

u = exp ^te + ^t^^ f{x + t^). 

By means of variables 

y = X + , y — hy = x — h~ + , 

y + h^ = x + h^ +t^, y + Ay = x + Ax + (t + r)^ 

we get the following system for the invariant solution of system (|7.19() : 



fiy + Ay)J exp^--Ay^-r(2y + Ay) + -r^ 



hy +hy y fiy) J hy y fiy) J J' 

where Ay can have one of the following values 

Ay = -hy, or (7.23) 

A solution of the above system with one of conditions (|7.23|) let us find the invariant 
solution for the operator Y4. 

5) Expressions , t~'^u, j and ^ are invariants of the operator I5. Let us search 
a solution of the difference model in the form 
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In variables 

X , _ X — h~ + X + Ax 

we get the following system of equations: 

/T—, ^ . > 2a / h- ( Hy + h+) \ ^h* ( f(y-h-) ' 



Here Ay can have one of the values determined by conditions H7.23|l and a is a constant 
from the condition a = j which determines an invariant time spacing. This condition can 
be found if we look for a time spacing r = g{t) which is invariant with respect to the 
considered operator Y5. 

6) For the operator Yq we have the following invariants: 

{4:f + l)^/'^iiexp + - arctan(2t) , 



4^2 + 1 , Ax 4*2 + 1 
+ At, 4t. 

T XT 

We look for a solution of the difference model in the form 
u = (tf + e.p (-^ - I „ct.n(2*)) / 

Involving new variables 

X X — h~ 

I X + h'^ X + Ax 

y + hZ = I , y + Ay = — , , 

" VWTT ^4(^ + ^)2 + 1' 

the system of equations (|7.19|) can be presented in the form: 

X exp ^carctan - b{y'^ + {y + Ayf) + 2\/62 + ly(y + Ay)J 



hy+hy \hj \ f{y) hy \ f{y) 
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where Ay has one of the values of (|7.2l-{j) . b is the constant from necessary condition of 
invariant grid existence 

4^2 + 1 
2b = U + —. 

T 

Therefore, the obtained reduced systems of equations determine the optimal system of 
invariant solutions for the difference model of the liner heat transfer equation. It means 
that each invariant solution can be found by transformation of a solution from the optimal 
system with the help of the corresponding element of the group. As we mentioned before 
the invariant time mesh for the new solution is obtained from the time mesh of the solution 
from the optimal system with the help of the same group transformation. For example, 
the transformation corresponding to the operator Xi with the value of the parameter —to 
gives shift in time t = t — tQ. Since 

Ad [exp(-to^i)]>5 = X, = Y5 + 2toXi, 

the action of this transformation transfers the invariant solution with respect to the op- 
erator Y5 into the solution which is invariant with respect to the operator X^,. By this 
transformation the spacing j = a is transformed into the spacing = a. 

Example of an exact solution. Among all group invariant solutions for the difference 
model (|7.19j) there is one interesting solution which can be integrated exactly |2|. This is 
the solution invariant with respect to the operator 

2toX2 + X3, to = const, (7.24) 

namely the solution 

considered on the mesh 

X. = x° (^) , (7.26) 

where are space mesh point at t = to- the case tg = we get the well known 
fundamental solution of the linear heat equation. Note that it has a "singular" mesh. 

Let us show how this solution can be obtained form the optimal system of the invariant 
solutions. From 



Ad [exp(eX5)]yi = X2 + 25X3 

we see that the solution (|7.25)) can be obtained from the solution invariant with respect 
to operator Yi by the transformation Ad [exp(eX5)] with e = If we take the original 
solution on the orthogonal mesh which is uniform in space and has the following special 
time spacing on the interval [0, to]: 

u{ = C, Xi = ih, i = 0,±1,±2, . . . , tj = i, j = 0,1,2,... 

to+JT 
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then the proposed transformation provides us the solution (|7.25|) on the uniform space 
mesh ()7.26() and uniform time mesh tj = jr. 

Thus, we see that difference model (|7.19|l inherits both the group admitted by the 
original differential equation and the ability to be integrated on a subgroup. 

4.5. The way to stop a moving mesh. The obtained difference models have adap- 
tive nonorthogonal grids. We can find a way to stop the moving mesh, i.e., an exchange of 
variables which orthogonalizes the mesh. The differentiation operator of Lagrange type ^ 
can be presented in the following form 

- = Dt- 2^D„ 
at u 



where 



^ d d 

ot au 



^ d d 
Dx = — + u^— + • 
ox ou 



The operator ^ in contrast to the operators Dt and does not commute with the 
operators of total differentiation with respect to t and x: 



UxUt 



U 



ul 



It is necessary to find an operator of total differentiation with respect to a new space 
variable s such that 



dt' 



0. 



(7.27) 



The last commuting property is possible if we involve a new dependent variable p > 
(density) ^Hl- The operator Dg = satisfies (|7.27j) if p holds the equation 



pt-2p[—-^\- 2-px = 0. 



u 



The new space variable s is introduced with the help of equations 

Ux 

St = 2p — , Sx = p. 



For convenience we can put initial data p(0, x) = 1. Then, s = x for t = 0. 

In the variables (t, s) the heat transfer equation will get a form of the system 



ut = P 



Uss - 2 — + ppsUs-, 
u 



pt = 2p- —-^ +2p^p,- 



u 

which can be rewritten in the form of conservation laws 



u 



P/t 



-2p^ 
u 



u 

P/t 



-pUs 



The space coordinate x is defined by the system of equations 

1 

p' 



xt 



Us 

u 



(7.28) 



(7.29) 



(7.30) 
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System (|7.28|) in the space of independent variables {t, s) and extended set of dependent 
variables {u, p, x) admits a group of point transformations determined by the following 
infinitesimal operators: 



d d ^ d d 

-^2 = -^■^ = 2t— xu-—, 

at ox ox Ou 



X4 = 2t^ + s-^, = 4^2 A + 4tx^ - {x^ + 2t)u-^ - 4tp^, 

ot ox OS ot ox Ou op 

= X-=/(4 + ,/'(.)|-, (7.31) 

where /(s) is an arbitrary function of s. 

In the independent variables (i, s) operators Xi-Xq are operators of Lie algebra fac- 
torized by the operator X* . Condition H2.6() of grid orthogonality and condition 1)2.5(1 of 
space grid uniformity hold and it gives an opportunity to construct a difference model 
which is invariant with respect to operators Xi~Xq on the orthogonal grid. 

Let us write system ()7.28p . ()7.3U() in the form of differential invariants. In the space of 
variables (t, x, s, u, p, dt, dx, ds, du, dp, Us , Ps,Xs,Uss) there are five invariants: 

Ug \ ids^'^ 
Jl = XsP, J2 = [dx + 2p—dt I , J3 — 




u ' ' p^dt ' 

Psds ^^3 f Uss u'^s\ o 2 _ 



U U^ J U 
\ f dx\^ ^ 2 f Uss \ . „ Us 



With the help of these invariants we rewrite system (|7.28() . I|7.3U() as 



ut = p^ 

P. = f ^ - 4 V 'ip'Ps-, = -2p^ (7.32) 











Uss 






U ) 






f Uss 




\ U 


U^ 



U u 



with the constraint equation Xs = - 



p 



Now we can find a system of equations which approximates (|7.32() and is invariant with 
respect to the set of operators (|7.31() . We can use the six-point invariant stencil (Fig. 4). 
We have the following invariants for the set of operators (|7.31|) where the operator X* is 
changed on its difference analog 

in the corresponding to the chosen stencil space {t,t, s, hf , , x,x, , , , ,u,u, 
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u+, u-,u+, u- , p, p, p+, p-,p+,p-): 



hx hx ^ hj u \A T 

1 /i+2 /In, 1 



"x ' "x 

'5 = T 



^In— + —In — , 

4 r h+ + h \ht u hx u J 

4 T ht + hx \ht u hx u J 

h~ , U: ht , 
^In— ^In — , 

hi u hx u J 

h~ ^ u+ ht ^ U-\ 
^In^ - ^In— , 
ht ^ hx u J 

T — I — T — ^ T — P+ T — ^t J 

-ffl — r -'in — — , -111 — , -'12 — ■ I ; -ilS 




P+ phi p-hx 

With the help of these invariants we can write difference model in the form of the 
following system of evolution difference equations (we present here only one invariant 
difference model which corresponds to the system ()7.19() in variables {t,x)): 



Ax = 2t- 



hs P In _i_ K P- l„ u- 

'h+P- u + h7 P u 

P P- 



n^^ ( lAx^\ ^ ^ ^In^ + ^ln^ 
- exp = 1 - 4r- ' 



uJ \ 2 T J hs I hs 

P P- 

ht 

i^x 

In the case of the uniform grid {ht = h~ = hs) this model could be simplified as follows: 

-pHn^+pl In^ 

hs{p + P-) 
ti\2 / lAx^X Arpp. 



u 



) '=''H"2— j = '-/5557:)('''"v+''-'-'ir)' 



ht 

The system ()7.28|) has only two dependent variables u and p and it can be approximated 
without involvement of the space variable x. However, Galilean symmetry and pro- 
jective symmetry are nonlocal in the coordinate system (t, s) and we need to consider 
the dependent variable x in order to have these symmetries. Constructing an invariant 
with respect to the set of operators H7.31|) difference model, we inevitably involve x into 
the difference equations. 

It is important to notice that in all cases moving meshes could be stopped by using 
the Lagrange type coordinate system (for an introduction of Lagrange type coordinate 
systems see, for example, 0^;). 
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5. If Q = Su, 5 = ±1, the equation 



ut = Uxx + 5u 



(7.33) 



can be transformed into equation (|7.12|1 by the change of variables 



u = ue 



(7.34) 



Reversing this transformation, one can get an invariant model for equation H7.33() from an 
invariant model for the heat transfer equation without a source. 
6. Q = 5 = const. The equation has the form: 



The case Q = 0, which was considered in point 4 of this section, is a partial case for 
Q = 5 = const, but the constant source can be excluded by the evident transformation 



It means that we can get a difference model for (|7.35|) from the model for (|7.12|) . 

8 Concluding remarks 

In the paper [10] the entire set of invariant schemes for ordinary difference equations of the 
second order was developed. There were shown that for some equations and symmetries it 
is necessary to involve nontrivial lattices, which are not uniform in a space of independent 
variable and should depend on solution. Such schemes are self adapted for any solution 
and they much exactly integrable, as its continuous counterpart. 

In the papers [2, 3, 8, 9, 11] several examples of invariant schemes for PDEs with 
two independent variables (KdV equation, nonlinear Schrodinger equation etc.) were con- 
structed. Again, there were shown that for some equations and symmetries it is necessary 
to involve nontrivial two-dimensional meshes, which are not uniform and rectangular in 
a space of two independent variables and should depend on solution. Such schemes are 
self-adapted for any solution, the meshes are evolutionary in time and these schemes have 
as many exact invariant solutions as their continuous counterparts. Moreover, for invari- 
ant variational cases invariant schemes have difference conservation laws as well as for 
continuous case. 

Thus, the applications of symmetry to difference equations led to the evolution of idea 
of possible meshes: from simple regular stationary meshes to self-adapted moving meshes, 
explicitly depending on solutions. Notice, that this idea evolution well corresponds to the 
big changes in numerical analysis, where the idea of self-adaptivity of schemes and meshes 
is in a broad fashion now. 

In this paper we developed the entire set of invariant difference schemes for the heat 
transfer equation 



Ut = u. 



'XX 



+ 5. 



(7.35) 



u = u — 5t. 



(7.36) 





for arbitrary coefficients K(u), Q{u) and for all special cases of the coefficients which 
extend the symmetry group admitted by equation (|8.1|) . 
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The main conclusion is that we have presented an algorithmic way to construct the 
invariant difference schemes (i.e. a difference equation and a mesh it is defined on) for all 
cases of underlying heat equation. 

Other conclusion is that symmetry preservation in difference schemes led to essential 
different discrete models. For different cases of coefficients K(u) and Q{u) taken in accor- 
dance with the group classification of equation (|8.1|) we have obtained different discrete 
models: for some cases K(u) and Q{u) we had to construct discrete models not for equa- 
tion (|8.1|) . but for the equivalent Lagrangian system: 



The consideration of the equation (|8.1() in the form of system ()8.2|) at the very beginning 
would provide us with the classification by functions if and ^. In that case there is 
one-to-one correspondence between the coordinate systems for continuous and discrete 
cases of the group classification. In particular, if the symmetry of the equation (|8.1|) does 
not require application of Lagrangian type moving mesh (evolving in time mesh) , then we 
have if = 0. In that case the classification of the system (|8.2|) and corresponding difference 
equations on orthogonal mesh can be carried out by means of 'ip{u,Ux,Uxx)- 
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